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ABSTRACT 

The zero-age main-sequence star AB Dor and the K dwarf component of the V471 Tau close binary 
have essentially identical rotation rates and spectral types. An analysis of their high resolution 
Chandra X-ray spectra reveals remarkably similar coronal characteristics in terms of both temperature 
structure and element abundances. Both stars show depletions of low FIP elements by factors of ~ 3, 
with higher FIP elements showing more mild depletions. No evidence for enhancements of very 
low FIP (< 7 eV) elements, such as Na, Al and Ca, as compared to other low FIP elements was 
found. The abundance anomaly pattern for AB Dor and V471 Tau is similar to, though less extreme 
than, the abundance anomalies exhibited by active RS-CVn-type binaries. While wc find statistically 
significant structure in the underlying differential emission measure distributions of these stars over 
narrow temperature intervals, this structure is strongly dependent on the lines used in the analysis and 
is probably spurious. Based on their X-ray similarities, we conclude that the exact evolutionary state 
of a star has little effect on coronal characteristics, and that the parameters that dominate coronal 
structure and composition are simply the rotation rate and spectral type. 

Subject headings: stars: abundances — stars: activity — stars: coronae — stars: late-type — Sun: 
corona — X-rays: stars 



1. INTRODUCTION 

The solar coronal abundance anomaly commonly 
known as the "FIP Effect" (or First Ionization Poten- 
tial Effect), in which low first ionization potential (FIP) 
elements (< 10 eV; e.g. Si, Fe, Mg) appear enhanced by 
average factors of 3-4, was already well established by 
the time the first observational clues to similar abun- 
dance anomalie s in stellar co ronae were uncovered in 
the 1990's (e.g. iFeldmanI 1199^ . These clues came from 
low resolution soft X-ray spectra {GINGA, ASGA, Bep- 
poSAX), together with moderate resolution extreme ul- 
traviolet (EUV) spectra obtained by the Extreme Ultra- 
violet Explorer (EUVE). The early stellar studies found 
evidence for a solar-like FIP Effect in some stars, but 
for the more active stars the abundances pointed toward 
metal-paucity rather than enhancement (e.g. .S. Drake, 
119961 and earlier references therein). 

In the last three years — the beginning of the Ghandra 
and XMM-Newton — era has seen the early hints of abun- 
dance anomalies fleshed out into an interesting array of 
diverse abundance patterns in which active stars appear 
to show signs not only of low FIP eleme nt depletion, 
but a lso of high FIP element enhancements l)Singh et al.l 
Il995|) . Studies of solar-like active sta rs confirin the s uspi- 
cions engendered by earlier work fe.g. l.T. Drakell99'^ and 
earlier references therein) that a transition from a metal- 
depleted t o a metal-rich cor ona occurs as the activity 
decreases ijGiidel et al.ll2002j) : this is now better char- 
acterised as a change f rom an apparently "in verse FIP 
Effect" to a FIP Effect. lAudard et aH (2003b') showed a 
similar transitions from an inverse FIP effect to the ab- 
sence of an obvious FIP effect for RS CVn binary stars 
with decreasing activity. Only one star studied in detail 
to date based on high resolution spectra appears to have 
a coronal composition demonstrably similar to that of its 



underlying p hotosphere — the subgiant Procyon (F5IV) 
l|Drake et a l.' 1995: ' R.aassen et al .' 20021. 

In the last decade, the com munity workin g on coronal 
abundances have argued (e.g. lDrakell2'003al and earlier 
references therein) that coronal abundances, when better 
understood, might provide new and powerful diagnostics 
of the physical processes underpinning stellar coronae. 
The emerging patterns of coronal abundance anomalies 
are telling us something about the dynamics structure 
and heating of coronal plasma; the challenge is in learn- 
ing to read these patterns. 

Aiming toward this goal, one question that arises is 
that of fine tuning of abundance patterns: how simi- 
lar are abundance anomalies i n stars that su perficially 
might have similar parameters? iDrakd l)2003b(l presented 
a differential analysis of the active binaries Algol and 
IIR1099 — two binary systems whose coronally-dominant 
components are of very similar spectral type and rotation 
period, but which have had quite different evolutionary 
histories — and found strong similarities in their quiescent 
X-ray spectra and coronal abundances. 

Here, we note that the young single active star AB Dor 
and the active component of the binary system V471 
Tau are both early K dwarfs and have rotation peri- 
ods of half a day. Superficially, we might also expect 
them to have very similar c oronal characteristics. How- 
ever, iStill fc HussainI 1)20031) have derived coronal abun- 
dances for V471 Tau based on low resolution ASGA CCD 
pulse height spectra which are significantly larger than 
those derived for AB Dor by other authors using simi- 
lar observational material, in addition to EUVE^^ XMM 



Newton and Ghandra spectra fRucinski et al.l '199 
iMewe et al._1 996; Ortolani et al._1998 ; Maggio et ai 
l200nt 



I et al. 2001; San z-Forcada et al. 2003). The 
apparent difference between the coronal abundances of 
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the two stars might be hinting that their coronae are in- 
deed fundamentaUy different. Such a dissimilarity might 
reflect the entirely different angular momentum evolution 
of the two stars: the rapid rotation of AB Dor is due to 
its relative youth, while that of V471 Tau is due to tidal 
spin-up from interaction with its close companion. 

In this paper, using Chandra High and Low En- 
ergy Transmission Grating spectrograph (HETGS and 
LETGS) observations, we present a comparative anal- 
ysis of the coronal X-ray spectra of AB Dor and V471 
Tau, with particular emphasis on their abundances. We 
first describe the two stars and briefly review earlier work 
then in fjSlwe report on the Chandra observations 
and data reduction. The methods used for a differential 
emission measure analysis together with results obtained 
are shown in 2] In f|Slwe discuss our results on the coro- 
nal and temperature structure and report our conclusions 
in Sg| 

2. PROGRAM STARS 
2.1. AB Dor 

AB Doradus (HD 36705) is a young, relatively bright 
(V=6.9) and rapidly rotating late-type star with spec- 
tral type K0-2V. It is an example of the very active 
cool stars that are just evolving onto the main se quence. 



Its age is esti mated to be 20-30 M yr (Guirado et al 
[ip97; Colhcr- Gameron fc Foind Il997t) . and llnnis et al 
l(l986.) suggested that AB Dor is member of the Pleiades 
based on its space velocity. The main ph ysical param- 
eters of AB Dor are shown in Table 'Ama do et al.l 
((200 1.) showed evidence of a solar-type activity cy- 
cle on AB Dor, with the starspot coverage increasing 
and decreasing with a period of 22-23 yr. Extensive 
Doppler imaging on this star has revealed the presence 
of a long term polar spot llCollier- Cameron fc UnruU 
1994[ .Donati fc C oUier C ameronI 1X997: Donati et aL, 
1999 1. The very fast rotation rate of AB Dor pro- 



duces strong magnetic fields which has been investi- 
gated thr ough Zeeman-Doppler Imaging by several au- 
thors fe.g. 'Donati fc Co llier Cameronll997tlDonati et al.l 
^99). Lim et al. (1994) presented measurements of the 
spectrum and polarization of the rotationally modulated 
radio emission on AB Dor. 

AB Dor was originally discovere d as a brigh t 
X-ray source, with Einstein, by iPakuU l)198lD . 
Since then it h as been exten s ively o bserved in 
the ultraviolet fRucinsk i et al.l 119951: iMoos et"al 



20001 iVilhu et all .200D and X-ravs llVilhu et al 
19871 ICollier Cameron etal] 1198^ IMewe et all 119961 



Kuerster et all 119971 iM aggio et all (20001 iGiidel et all 
20011) . iKuerster et al . (1997) found variability on 
short timescales but no long-term trend or cyclical 
variability of the X-ray emission. Their study showed 
that about 15% of the X-ray emission is rotationally 
modulated, indicating that at least part of the X-ray 
emission must have scale si z es below a ste l lar ra dius . 
Collier Cam eron ct al.l l)1988|) . iSchmitt et al] 1)19981 ) and 



Maggio et al... (,2000i) have reported X-ray flares on AB 



Dor. 

The photospheric abundances of AB Dor are rela- 
tively well studied; Vilhu ct al. l)1987t) found a solar- 
like metallicity, [Fe/II]=0.1±0.2, based on a high reso- 
lution spectroscopic study. The Fe photospheric abun- 
dance reported for the Pleiades by several authors (e.g.; 



IKing et all I2OO0I: IBoesgaardI IT989I: ICavrel et all IT981 
show a slightly super-solar value, 0.02<[Fe/H]<0.15. 
Similar values have been re porte d for other elements 
([O/H]==0.14 (Schu ler et all l200l . [Mg/H]=0.05 and 
[Si/H]=0.07 (King ct al.' "2000)). A low coron a l iron 
abundance, [Fe/H]=-0.95, found by IMewe et al J (CSal) 
using simultaneous observations w i th EU VE and ASCA, 
was confirmed by IGiidel et all \2QQl\i using XMM- 
Newton data and what appears to be an "inverse FIP" 
effect — that elements with high First Ionization Poten- 
tials (FIPs) (e.g., C, N, O, Ne, Ar) appear enhanced 
relative to low FIP elements (e.g., Mg, al. Si, Fe), oppo- 
site w hat is observed in the solar corona (e.g. iFeldmanl 
Il992|) . These values were 3-5 times lower than the 
photospheric abundance obtained bvlVil hu et al.l (|1987^. 
Sanz-Forc ada et al. (2003) reported coronal abundances 
for AB Dor based on XMM-Newton and Chandra spectra 
showing a similar Fe depletion. 

2.2. V47I Tau 

V471 Tau {BD+16 516) is a rapidly rotating rela- 
tively fai nt (V=9.5) detached ecl ipsing binary first ob- 
served bv lNelson fc Youn g^ ('197QD. The system is com- 
posed of a co ol main-se quence chr omospherically ac- 
tive K2V star (jCuinan fc Si on 1984) and a degenerate 
hot white dwarf (DAI. 5) separated by only 3.4 stel- 
lar radii. During the primary eclipse, the hot white 
dwarf acts as a beamiiig pro be for the K2 star's at- 
mosphere (jCuinan et alj[l98^ . The system is thought 
to belong to our nearest open cluster the Hyades based 
on Hipparcos proper-m otion and paralla x measurements 
(jProvencal et alJll996l Tde Bruiin e et all UoOll although 
this point has be en controversial (see discussion by 
iMartin et aDll997D . The basic properties of V471 Tau 
are summarized in Tabled 

As discussed by Paczvnskj l)1976ft . the V471 Tau sys- 
tem has probably been through a common-envelope (CE) 
stage, and the system is evolving to become a cataclysmic 
variable in which the K dwarf overfills its Roche lobe, 
leading to mass transfer onto the white dwarf. V471 
Tau is often considered to b e the prototypical " pre- 
cataclysmic" binary (see |Bond 198^; it is one of two 
known post-CE, pre-CV systems in the Hyades clus- 
ter (the other being HZ 9). Based on the same Chan- 
dra LETGS spectrum that we analyse in this paper, 
[Drake fc Sarna (2003) provided the first direct observa- 
tional evidence of the CE phase of V471 Tau, demon- 
strating that the coronal C abundance was depleted rel- 
ative to N as would be expected from the surface contam- 
ination of the K dwarf by the evolved and expanded enve- 
lope of the primary star. V471 Tau is currently the only 
known precataclysmic v ariable detected at radio wave- 
lengths l)Lim et allligg^ . 

Several observers have reported a strong, complex, 
and variable Ha feature that goes from absorption 
to emission as the system rotates (Young et al. 198^ 
iSkillman fc Pattersonl[T98^ . Other chromospheric lines 
such as the Ca II H & K and Ca II infrared triplet 
show similar behavior confirming that the activity is 
coming from the chromosphere of the K dwarf. Dif- 
ferent authors have presented evidence for magneti- 
cally induced activity on the K dwarf based on op- 
tical and photometric f l ares jYoung et al.l 119911) . X- 
ray flares ijYoung et al.1 11983(1 and radio flaring ac- 
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tivity llPatterson et all 119931: iNicholls k StorevI Il999|) . 
I.Tensen et alJ |)198(tD and lsion et a].l 1119981^ found a pe- 
riodic 9.25 minute modulation in the soft X-ray and 
EUV bands of V471 Tau, subsequent ly found in the 
optical band by Robin son et al.l 1)198 8'^ . It is now be- 
lieved that the variability is caused by rotational modu- 
lation of the magnet ic white dwa r f llClemens et alJll992l: 
iBarstow et al]ll993j) . iWheatlevI found that the 

rotation rate and X-ray luminosity of the K star in V471 
Tau are typical of K stars found in the Pleiades. 

The photospheric abundances for the Hyades open 
cluster, to which V471 Tau is thought to belong, are 
generally better determin ed than thos e for m ost coro- 
nally active field stars. iCavrel et alJ (|1985D reported 
an analysis of high signal-to-noise ratio Coude spectra 
for 12 Hyades dwarfs, together with solar spectra for 
comparison, in order to obtain accurate abundances for 
iron and other metals in the Hyades. They determined 
mean Fe and Si abundances slightly higher than the so- 
lar values, [Fc/Hl=0.12±0.03 and [Si/H]=-t-0.16±0.04. 
IVarenne fc Monier (1990,) derived abund ances for O I, 
[Q/H| =-t-0. 07 for a similar sam ple. Boes gaard fc Budgg 
l)1989tl and lVarenne fc Monied l|l999i) found, based on a 
sample of F stars in the Hya des, similar values for the 
[Fe/H]. iSmith fc Ruc3 l)1997|) reported a global metal- 
licity for the Hyades cluster of [Fe/H]=-t-0.13 based on 
values determined for i ndividual dwa rfs and giants from 
high resolution spectra. iSmithI l)1999() derived similar val- 
ues for the Hyades KO giants 7 Tau an d e Tau. Based 
on low resolution X-ray CCD spectra, iStill fc HussainI 
(120031) reported coronal abundances showing only a very 
mild Fe abundance depletion relative to photospheric es- 
timates, together with evidence for an inverse FIP effect. 
In their analysis Fe coronal abundance shows a factor of 
^2 depletion while Ne shows a factor of '^3 enhancement, 
both relative to photospheric values. 

3. OBSERVATIONS 

The Chandra HETGS observation of AB Dor was car- 
ried out in 1999 September using the AXAF CCD Imag- 
ing Spectrometer (ACIS-S), while the Chandra LETGS 
observation of V471 Tau was carried out in 2002 Jan- 
uary using the High Resolution Camera spectroscopic 
(HRC-S) detector. Both observations employed their re- 
spective detectors in their standard instrument configu- 
rations. The observations are summarized in Table [3 

Pipeline-processed (CXC software version 6.3.1) pho- 
ton event lists were reduced using the CIAO software 
package version 2.2, and were analyz ed using the IDL^- 
based PINTofALE software suite ijKashvao fc Drak^ 
EOOO). In the case of V471 Tau, processing included fil- 
tering of events based on observed signal pulse-heights 
to reduce background (Kashyap fc Drake 2000, Wargelin 
et al., in preparation). The analysis we have performed 
consisted of line identification and fitting, reconstruction 
of the plasma emission measure distribution including 
allowance for blending of the diagnostic lines used, and 
finally, determination of the element abundances. 

The white dwarf component of V471 Tau exhibits 
a strong thermal continuum at long wavelengths (> 
50 A), which is clearly seen in the LETG spectrum. 
This was studied in previous EUV observations (e.g. 

^ Interactive Data Language, Research Systems Inc. 



IBarstow et al.lir992() . We note, however, that this con- 
tinuum contributes no significant flux shortward of 50A. 

Fig. n shows the Chandra X-ray spectra of AB Dor 
(black shade; HETG-t-ACIS-S) and V471 Tau (grey 
shade; LETG-I-HRC-S) in the limited wavelength range 
3-27 A which contains the prominent lines of Na, N, O, 
Ne, Mg, Fe, S and Si. The strongest coronal lines are 
identified. The thick black line shows the AB Dor spec- 
trum after smoothing in order to match the resolution 
of the V471 Tau spectrum. The spectra of both stars 
show a remarkable similarity of lines, both in terms of 
which lines are prominent, from H- and He-like ions and 
the broad range of charge states of Fe, and in their in- 
tensities. Figure^shows that at shorter wavelengths the 
spectrum is dominated by emission lines that we inter- 
pret as coronal emission from the K dwarf compo nent. 

According to the ephemeris of Guin an fc RibasI (|2P01), 
the LETGS observation of V471 Tau started at phase 
4) = 0.05 and ended at = 1.97, where (p — when the 
K dwarf component is nearest the observer. 

4. ANALYSIS 

4.1. Photometry 

Before commencing spectral analysis, we first checked 
for flare activity that could affect not only the shape 
of the differential emission measures (DEMs) but might 
also be accompanied by detectable changes in the coro- 
nal abundances of the pl asma that dominates disk- 
averaged spectra (see. e.g. lOttmann fc SchmittI 119961: 
i Mewe et al.ll997tlOrtolani et al.ll998l: | Favata fc Schmitt 
'1999>'Favata et a,].l200nHMa,ggio et a,].l200nHGiide] et al. 
2001). Light curves for the AB Dor and V471 Tau obser- 
vations were derived using the 0th order events extracted 
in a circular aperture, in the case of V471 Tau, and using 
the dispersed photon events and excluding 0th order (to 
avoid photon pile- up) in the case of AB Dor. Events were 
then binned at 100s and 200s intervals for AB Dor and 
V471 Tau, respectively. Fig. [21 shows the light curves for 
both observations; note that both light curves are rela- 
tively flat and devoid of large flare activity. 

Several authors have reported coronal r otational 
modulation o n active stars (e.g. iDrake et al.. 199^ 
lAudard et aLir20 01: Garcia- Alv arez et all 12003(1 . includ- 
ing AB Dor jKuerster et al.«il997j) . However, as is shown 
in Fig. 121 (bottom panel), no obvious modulation of the 
coronal emission with orbital phase was seen on either 
star, though in the case of AB Dor, the observation only 
covers slightly more than one rotational period. In the 
case of V471 Tau, modulation due to modification of the 
underlying chromospheric structure might be expected, 
for example, if the part of the K dwarf corona irradiated 
by the white dwarf were affected by the UV radiation 
field of the latter. This is apparently not the case to any 
significant extent. A similar conclusion has also been 
reached bv IWheatlevI (|l998). 

We conclude that both AB Dor and V471 Tau obser- 
vations are representative of the stars during times of 
quiescence, and therefore treat the observations in their 
entirety for the remainder of the analysis. 

4.2. Spectroscopy 

Spectral line fluxes for both AB Dor and V471 
Tau were measured by fitting modified Lorentzian ("/?- 
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profile") functions of the form 
F(A)=a/(l + ( 



A- Aq 

r 



Y) 



(1) 



where a is the amplitude and T a characteristic line 
width. For a value oi (3 — 2.4, it has been found 
that this function represents the Chandra transmission 
grating instrumental profile to within photon co unting 
statist ics for lines with a few 1000 counts or less ijDrakd 
l2003c^) . Spectra line fitting was undertaken using the 
FITLIN ES utility in the PINTo fALE spectral analysis 
package ijKashvap fc Drak3l200nD .^ In the case of ohvi- 
ous blends, as in the Ne IX-Fe XIX 13 A region, we have 
performed multi-component fitting. 

While strong lines can be easily identified, the sum of 
weak lines, each of which can be undetectable, can pro- 
duce a "pseudo-continuum" . The true continuum level 
was set using the spectral regions 2.4-3.4, 5.3-6.3 and 
19.0-20.0 A, which we judged to be essentially "line-free" , 
based on both visual inspection and the examination of 
radiative loss models. These continuum points were then 
used to normalize model continua computed using a test 
Differential Emission Measure taken from a n analysis of 
a typi cal active star (AU Mic in this case, iDrake et al.l 
l)2004l) V In principle, the fiuxes should be re-measured 
once the final DEM has been determined; in practise, 
there was little difference in the shapes of trial and final 
continuum models so that this additional iterative step 
was not required. 

The emission measure distribution and abundance 
analy sis employed the CHIANTI datab ase (iDere et al 
120011) and the ionization balance of iMazzotta et al 
()1998|) . as implemented in the PINTofALE software 
package (Kashvao & Drake' 12000(1 . Table |3| shows the 
measured fluxes of the emission lines identified and used 
in this analysis. We have computed correction factors 
(CF) for line blends by computing complimentary error 
functions for each contaminant (i.e. Gaussian profiles 
are assumed). A set of interesting lines and a set of 
contaminant lines, which contain necessary line informa- 
tion (emissivity, wvl, z, etc.), are used. The procedure 
identifies which lines are possible contaminants for each 
analysed line. The contaminants must lie within two 
sigma on either side of the line and be greater than a 
user-defined threshold. The CF are defined as follows: 

/(A.) + E,^./(A.) 

where I(Ai) = intensity of the line being contaminated 
by lines j ^ i- CF for line blends for the studied lines 
range from 0.89 to 0.99. Fig. |21 shows the comparison 
of our observed and blend corrected fiuxes for both AB 
Dor and V471 Tau. Note that the fluxes are corrected 
for distance and exposure time. 

4.3. Differential Emission Measure 

We adopted the following formalism for the Differen- 
tial Emission Measure. The line flux, for optically-thin 
plasma, for the transition of an ion between levels i j 
can be written as: 

F,, = AK,, [ G,y (ne,T)$(r) dT erg cm'^ s^^ (3) 

^ PINTofALE is freely available from 

|http://hea-www.harvard.edu/PINTofALE/| 



where A is the abundance, Kji is a known constant which 
includes the frequency of the transition and the stellar 
distance, ne is the electron density, G^ is the contribu- 
tion function and T is the electron temperature. The 
kernel <I>(r) is commonly known as the differential emis- 
sion measure (DEM) and is a measure of the amount of 
emitting power — correlated with the amount of emitting 
material — as a function of temperature in the coronal 
plasma. It is formally defined as: 



dV{T) 
dT 



(4) 



Two different approaches are commonly employed in the 
derivation of the DEM in the corona, namely, line-based 
methods and global fitting techniques. For the latter the 
parameters are determined by a minimization of the 
observed and trial model spectra. We prefer the former 
method for the greater control over the spectroscopic di- 
agnostics entering the analysis that it allows. 

A given line flux depends on both temperature and 
on the abundance of the clement in question. For this 
reason, methods wishing to take advantage of bright 
lines from a number of different elements generally have 
to perform simultane ous modelling of the DEM and of 
abundances (see, e.g. iMaggio et alJl2004l: lArgiroffi et alJ 
^04). However, the ratio of two emission lines from 
ions of the same element is independent of the abun- 
dance of the chosen element. We have therefore devised 
a method that uses line ratios instead of line fluxes di- 
rectly. A similar method was developed independently 
by Schmitt & Ness (2004) who also discuss the benefits 
of such an approach. Our basic set of diagnostics com- 
prises the H-like/He-like resonance line flux ratios for 
the elements O, Ne, Mg, and Si, line ratios involving 
Fe XVII, Fe XVIII and Fe XXI resonance lines, and mea- 
surements of the continuum flux at points in the spec- 
trum that are essentially free of lines (the same regions 
noted above with regard to spectral line intensity mea- 
surement). Our m ethod is described in more detail by 
l|Drake et al.ll200^ . 

The motivation underlying the development of this 
method lies in a desire to survey coronal abundances and 
trends in temperature st ructure independen t from issues 
regarding atomic data. IKashvap & Drak3 l(1998tl have 
demonstrated that a reconstructed DEM can be quite 
sensitive to the lines used for the reconstruction, and 
also on the atomic data adopted for those lines. Our ba- 
sic line set comprises the brightest lines in stellar coronal 
spectra and are easily measured in essentially all well- 
exposed Chandra grating observations of stellar coronae, 
such that star-to-star variations in diagnostic lines used 
can be avoided. The H-like and He- like ions are also likely 
to be represented by more accurate theoretical line inten- 
sities than lines from more complex ions. Similarly, the 
theoretical description of free-free and bound-free con- 
tinua used in our DEM analysis is thought to be well- 
understood. In particular, we note that the continua in 
these hot coronal sources tend to be dominated by free- 
free emission. 

4.4. Algorithm and Line Selection 

In order to obtain the differential emission measure, 
^(r), we have performed Markov-Chain Monte-Carlo 
analysis using a Metropolis algorithm (MCMC[M]) on 
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the s et of suppUed Une flux ratios ijKashvap fc Drak3 
Il998j) . Advantages of this method include the abil- 
ity to estimate uncertainties on the derived function 
$(T), and the avoidance of unnecessary smoothing con- 
straints. Eq. 01 defines the DEM as a continuous func- 
tion of temperature T. The MCMC[M] method yields 
the emission measure distribution over a pre-selected 
temperature grid, where the DEM is described by a 
simple histogram. In our case, a set of temperatures 
T„, with A log T[K]=0.\ and ranging from log r=6.1 
to log T[K]^%.0, define the $(T„). The derived $(T„) 
is only reliable over a certain temperature range if we 
have enough lines with contribution function G{Tmax)^ 
G(T„)._ 

The inversion problem is ill-conditioned which can- 
not be overcome by any fit procedure criteria applied 
ijCraig fc Browrll976|) . Nevertheless, in order to "avoid" 
the problem of mathematical ill-conditioning, DEM re- 
construction techniques often rely on modeling the DEM 
as smooth splines or polynomials, and are thus artifi- 
cially constrained to be smooth. We do not impose ar- 
bitrary smoothing restrictions, but instead let the algo- 
rithm find the best possible m odel allowed by the data 
(e.g.; iKashvap fc Drakj 119981) . while applying a vari- 
able local smoothing to the DEM. This local smooth- 
ing is based on the diagnostic information available: if 
there are many spectral lines in a given temperature 
range with different temperature sensitivities, smooth- 
ing is essentially switched off; instead, more smoothing 
is applied where diagnostic information is more sparse. 
This technique often allows solutions containing struc- 
ture over small temperature ranges, whose formal signif- 
icance may be judged by the sizes of the error bars that 
the algorithm also produces. While the reality of strongly 
peaked features in stellar DEMs is still in question, other 
techniques have found evidence for sharply peaked struc- 
tures in some coronae, as illustrated, for example, by the 
emission measure d erived for Capella using EUVE data 
(iDupree et alJll993j) and for the Sun using Yohkoh data 
l)Peres et al.ll200 ll)7 Although we have found a "good" 
solution for the DEMs, by using the local smoothing cri- 
teria, one has to keep in mind that we cannot a priori 
exclude other solutions in the parameter space that de- 
scribe the spectra equally well. 

Table O shows the lines used for the DEM(T„) recon- 
struction. In the X-ray range analyzed, the lines with the 
coolest peak formation temperature are the resonance 
hues O VII and N VII, log T[K]=6.3, while the hottest 
peak formation temperature is the resonance line S XVI, 
log T[K]=7A. Based on the lines we use in our analy- 
sis we are able to obtained a well-constrained DEM(T„) 
between log T[K]=6.3 and log T[K]=7A; larger uncer- 
tainties are obtained outside that range. 

In order to probe the consistency of the results, de- 
pending on the line ratios used to derive the DEMs, we 
chose four different sets of H-like, He-like and highly ion- 
ized Fe lines. Set I contains H-like and He-like lines only 
(O, Mg, Ne, Si), Set 2 contains O, Mg, Ne, Si and S, 
Set 3 contains O, Si, S, Fe XVII, Fe XVIII and Fe XXI 
and Set 4 contains O, Ne, Mg, Si, Fe XVII, Fe XVIII 
and Fe XXI. Set I repres ents the bas i c core of H-like 
and He-like lines used by ' Drake et al.l l)2004() to derive 
DEMs. Set 2 adds sulfur to the previous set of lines, 
which allows us to increase the upper temperature range 



up to log T[K]=7.4:, although the significance and signal- 
to-noise for these lines in our spectra are low. Set 3 con- 
tains H-like and He- like lines that constrain the structure 
of the DEM in the upper and lower temperature range, 
while the iron ions give information about the structure 
of the DEM in the intermediate temperature range. Set 
4 contains all the ions analyzed. All the sets contain 
O VIII, which constrains the DEM in the lower tempera- 
ture range, log T[K]=6.3. Note that the electron density 
(ug) adopted for Eq.^will not produce significant differ- 
ences in the derived DEMs since the chosen lines are not 
Ue-sensitive. 

The MCMC[M] code returns the DEM, among any 
user-defined number of possible DEMs that are gener- 
ated based on the MCMC method, that best reproduces 
the observed fluxes. As noted above, one product of 
our MCMC[M] method is the statistical uncertainties in 
$(T) evaluated at the 68% confidence level. One has to 
bear in mind that the error bars for consecutive bins in 
DEMs are not independent, owing to the inevitable cross- 
talk that arises as a result of the line contribution func- 
tions themselves stretching over several bins. While the 
DEM method imposes a degree of local smoothing, dic- 
tated by the number of diagnostic li nes with significant 
emission in each temperature bin fsee IKashvap fc Drakd 
1998, for details), this smoothing is not always suffi- 
cient to attenuate "high frequency oscillations" in the 
solution — the artifact whereby overly high values of the 
DEM in a given temperature bin can be compensated for 
by a commensurately lower value in adjacent bins. 

Our final DEMs are corrected fo r the distances of the 
stars: we have adopted £)=14.9pc iGuirado et al.lll997[l 
and 15=48 pc l|de Bruiine et alJ 120011) for AB Dor and 
V47I Tau respectively. 

4.5. Element Abundances 

Once the DEM, $(T), has been established, we can 
evaluate the abundances of any elements for which we 
have lines with measured fluxes. The abundance is then 
a simple scaling factor in the integral equation [SJ and is 
given by the value this factor has to assume in order to 
match the measured flux. In this way, we derived values 
for the coronal abundances of the elements O, Ne, Mg, 
Si, S and Fe, though only upper limits were obtained 
for the abundances of S for V47I Tau owing to the lack 
of signal in both Hc-like and H-like features. According 
to previous works the photospheric abundances for AB 
Dor and V47I Tau are similar to the solar ones (se e i )2.I| 
and i )2.2|l . Therefore, we will use the|G revesse fc S au val 
((T998ijl solar photospheric composition but will refer to it 
as stellar photospheric values. Note that we have used 
the new solar p hotospheric oxygen abunda nce, 0/H = 
8.69, derived bv lAllende Prieto et in (l200l . 

We also used th e temperature-insen sitive abundance 
ratio diagnostics of iDrake et al.l l)2004l) as an additional 
check on our values obtained using the DEM. These 
are ratios of formed by combining two sets of lines of 
two different elements, constructed such that the com- 
bined emissivity curves of each set have essentially the 
same temperature dependence. The resulting ratio of 
measured line fluxes then yields directly the ratio of 
the abundances of the relevant elements, independent 
from the atmospheric temperature structure. A similar, 
though less sophisticated, method has been previously 
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used in solar coronal observa tions (see, e.g. iFeldmaiJ 
119921 iFeldinan fc LaminJ20f)nD and in stellar EUVE ob- 
servations ijPrake et alJll995l I1997D . 

5. RESULTS AND DISCUSSION 

In order to verify the propriety of our DEM and abun- 
dance techniques, we have compared observed and mod- 
elled line fluxes vs ionic species in Fig2|(upper panel) and 
vs Tmax (lower panel) for the four line sets and for the 
final DEMs and abundances. Open symbols and filled 
symbols show the values for AB Dor and V471 Tau re- 
spectively. The great majority of the predicted line fluxes 
based on our final models are within 10 % or so of the 
observed values, with only four lines being dis crep ant by 
more than 20% . We address these further in iiS.ll below. 

Fig 13 and Fig |H1 show the observed and synthetic 
spectrum of AB Dor and V471 Tau. The synthetic 
spectra were computed using line emissivities from the 
CHIANTI database, and were obtained with the de- 
rived DEMs and abundances, for the four line sets. We 
also computed synthe tic spectra using e missivities from 
the APED database l|Smith et al J 120011) . and the DEM 
and abundances derived from line set 1 for compari- 
son. For these simulations, we have use an electron den- 
sity of ne=5 lO^^cTO"^ for AB Dor and ne=2 10^''cm~'^ 
for V471 Tau (Testa et al. 2004), with an interstellar 
medium column density of n^f =1 10^'^cm^^. The exact 
density and very low ISM column makes no practical dif- 
ference to the synthetic spectra. All the spectra are re- 
markably similar with all the synthetic spectra showing 
good qualitative agreement with observations for both 
stars. 

5.1. Temperature Structure 

The reconstructed DEMs for AB Dor and V471 Tau, 
together with estimated la uncertainties, are shown in 
Fig. 13 The DEMs are all quite similar in overall shape 
and normalisation for the four sets of lines used. A peak 
at log T[K]^7.0±0.1 is observed in all the reconstruc- 
tions independently of the set of lines used. The DEMs 
for AB Dor also exhibit secondary peaks at logT[K]~ 
6.6±0.1 and logT[K]- 7.4±0.1. The DEMs obtained by 
using the line sets with no iron ions (sets 1 and 2), show 
sharper increase and decrease for the lower and higher 
temperatures compared with ones obtained by including 
iron lines (sets 3 and 4), which show more prominently 
defined peaks. Of these, it is interesting that the DEMs 
for both AB Dor and V471 Tau corresponding to line set 
4, containing all the ions, require the most structure. 

The DEMs for both stars show little evidence for sub- 
stantial amount of plasma at temperatures higher than 
log r[K]=7.4. Nevertheless, the corona of AB Dor seems 
to be, based on the DEMs, slightly more active than that 
of V471 Tau. The slopes in the derived DEMs to the 
low temperature side of their maxima follow very ap- 
proximately the relation EMocT^, which is steeper than 
the EMocT'^/^ suggested for coronal lo ops with constant 
cross-section and uniform heatin g fe.g.: ICraig et al.ll978t 
IJordanllT980t IPeres et al.ll200(]|) . If the sharp peaks ob- 
served in the DEMs are real, it would imply that the 
corona must contain essentially isothermal structures 
(see also lDrake et a l. 2000:' Scelsi et al.ll2004D . 

Our results for AB Dor are comparab l e with those 
reported recently by iSanz-Forcada et al] l)2003f) based 



on yiMM-Newton and Chandra observations. They ob- 
tained a three-peak DEM with a more smooth increase 
for the lower temperatures, log T[K]<6.6, and sharp de- 
crease for the higher tein peratures, log T[K]>7.4. The 
ISanz-Forcada et al.l lj2003 t) line-based DEM also shows a 
maximum at log T[K]~7.0±0.1. To our knowledge, there 
exist no other DEM analysis of V471 Tau by other au- 
thors with which to compare our temperature structure 
results. 

How reliable are the structural features observed in 
the DEMs? Kashvap fc Drake (1998) demonstrated that 
solar quiet and active region DEMs are not stable to 
changes in the lines used in the analysis. This occurs 
because of a breakdown in the assumptions of optically 
thin plasma analysis, or (or in addition to) because of 
uncertainties in the underlying atomic data that are not 
included in the analysis. We note in passing that mea- 
surements of the O VIII and Ne X Lya/Ly/3 ratios for 
AB Dor are in perfect agreement with the optically thin 
theoretical values (Testa et al., in preparation), so that 
we do not expect photon loss through, e.g., resonance 
scattering to be important. 

The principal lines used in our analysis are simple H- 
like and He-like transitions for which atomic data should 
be re liable. Th e exceptions in line sets 3 and 4 are the Fe 
ions. IDoron fc' Bchar (2002) and Gu (2003) reported new 
Fe L-shell calculations indicating that rate coefficients 
for all n— >2 (3<n<5) transitions may differ from earlier 
calculations up to 50%. Moreover, these authors em- 
phasised the importance of including in the level popu- 
lations the indirect processes of radiative recombination, 
dielectronic recombination, and resonance excitation in- 
volving neighbouring charge states. It seems likely that 
problems associated either with the Fe line emissivities 
themselves, or with the underlying ionization balance, 
might be to blame for the differences in the temperature 
structure obtained from the different line sets. Regarding 
the ionization balance, we also note that dielectronic re- 
combination (DR) tends to dominate the recombination 
rates for the Fe ions near their temperatures of peak pop- 
ulation. In contrast, for H-like ions, DR is not an issue 
and one seems justified in having more faith in the pop- 
ulations of these ion s. The case of He-l ike ions has been 
discussed recently bv lSmith et alJ l|200H) based on the ex- 
ample of O VII. Unlike the excitation of the (ls.2p) 3Si 
and (ls.2p) Pi_2 levels that form the upper levels of the 
prominent forbidden and intercombination lines, respec- 
tively, and which tend to be dominated by dielectronic 
recombination and cascades, the (ls.2p) IPi upper level 
of the resonance line is excited almost entirely by colli- 
sions. Thus, the He-like resonance lines should also be 
quite accurately represented in recent theoretical models 
and should provide reliable diagnostics. 

As we noted earlier in a small fraction of our pre- 
dicted line fluxes based on final DEMs and abundances 
differ to those observed by 20 % or slightly more (Fig^J. 
Conspicuous among these are the O VII lines in AB Dor. 
This small discrepancy is likely due to the lack of strong 
temperature constraints for the cooler temperatures — a 
similar, though smaller, effect is discernible at the high 
temperature end from S XVI. One other conspicuous 
point in Fig 0] belongs to Ne IX in AB Dor for line set 
4 that includes Fe lines in addition to the H-like and 
He-like ions. Since this No IX predicted flux is the only 
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discrepant point among the three line lists in which it 
appears, we conclude that there is no DEM solution that 
can accommodate both the Fe and other lines simultane- 
ously to within their statistical uncertainties. In essence, 
the Fe XVII and Fe XVIII lines, that are formed at tem- 
peratures similar to those of Ne IX and Ne X, are pulling 
the DEM to a solution that is unfavourable to the Ne IX 
line. 

We have argued above that atomic data uncertain- 
ties are a likely culprit for this situation. However, an- 
other possibility is that the coronal abundances are non- 
uniform as a function of temperature. This would also be 
seen in the differences between DEMs derived by using 
set 1 and set 2, in which we use Mg and Ne ratios for the 
low temperature range, and set 3 in which we swapped 
those ratios with Fe ones. The latter results in a much 
lower DEM at logT ^ 6.6, near where Fe XVII ions 
reach maximum population. To retain a similar DEM to 
those obtained from line sets 1 and 2 at this tempera- 
ture would require an increase in the Fe abundance at 
this temperature relative to the value at higher tempera- 
tures. Equally, similar agreement of the DEM's from line 
sets 1, 2 and 3 might be achieved by a decrease in the 
Ne abundance toward higher temperatures. These argu- 
ments are qualitative and we have insufficient spectro- 
scopic diagnostic information to investigate abundance 
variations with temperature in a rigorous fashion. We 
note, however, that such abundance variations might not 
be unreaso nable in the Hght of a recent survey of coronal 
densities bv lTesta et al.l ("2004") that clearly demonstrates 
that the cooler T ~ 10*^ K corona resides in quite differ- 
ent structures to the hotter, T ~ 10^ K plasma that 
appears t o be universally ch aracterised by much higher 
pressures. IFeldman &: Lamm g (2000) found that for the 
solar corona the FIP effect is larger for T ~ 10^-^ K than 
at T ~ lO^ is K. 

Regardless of the true cause of the differences be- 
tween the DEMs obtained from different line sets, we 
can conclude that small-scale structure obtained in coro- 
nal DEMs is not reliable: such structure sensitive to 
the particular diagnostics chosen and does not neces- 
sarily reflect the true underlying temperature structure. 
This conclusion echoes the statements made earlier by 
IKashvap fc D?a kc (19980. 

The DEMs of V471 Tau are more similar among the 
different line sets than those of AB Dor. However, de- 
spite of their quite different evolutionary histories, we 
conclude that their coronae show a strong similarity in 
temperature structure. 

5.2. Coronal Abundances 

Our results for the clement abundances of AB Dor and 
V471 Tau are summarized in Tabled] We only list sta- 
tistical uncertainties here. There are of course other un- 
certainties in the atomic data used for the analysis, and 
in the calibration of the Chandra instruments. Atomic 
data uncertainties are probably of order 30% for Fe lines 
(e.g.. IDrake et alJll995^) . but are probably slightly less 
than this for H-like and He-like lines. The final uncer- 
tainties in our derived abundances are difficult to assess 
rigorously, but are likely of order 0.1 dex. 

The abundances derived are quite consistent for the 
different line sets considered: the abundances for each el- 
ement agree even within the statistical errors except for 



oxygen in set 3 of AB Dor. We also notice that the abun- 
dances for oxygen and iron for both AB Dor and V471 
Tau derived from set 3 are slightly higher than those de- 
rived from the other sets, an effect due to differences in 
the temperature structures obtained (see ti5.1|) . 

Table[Slshows the abundance ratio, from He-like and H- 
like ions of light elements, that are relatively insensitive 
to temperature (see t|4.5|l . Those values are in agreement 
with the ones derived from the DEMs (Table^J, with the 
exception of the ratios containing iron. Atomic physics 
uncertainties in the iron ions might affect the values ob- 
tained for the ratios [Ne/Fe] and [Mg/Fe]. 

In Fig.|Slwe have plotted the derived abundances, rel- 
ative to the stellar photospheric values, in order of the 
element FIP. AB Dor and V471 Tau clearly share similar- 
ities in coronal composition: both show the same general 
pattern characterised by a depletion of low FIP elements 
by a factor of '--^ 3, evidence for slightly smaller depletions 
of intermediate FIP elements S and O, and a clear en- 
hancement of Ne relative to elements of lower FIP. This 
pattern is similar to, though less extreme than, those seen 
in other ac tive stars and RS CVn-type binari es (see, e.g., 
reviews bv lDrakell2003aHAudard et alJl2003b|) . That the 
abundances found here represent real metal depletions 
in the coronae of both AB Dor and V471 Tau is beyond 
doubt. As discussed in ^ both stars have solar-like, or 
slightly super-solar, metal abundances and this is irrec- 
oncilable with the much lower coronal abundances of low 
FIP elements. 

In Table |S1 we list the abundance results from earlier 
studies of AB Dor and V471 Tau (Mcwc et al. 1993 
l^tolani et al. 1998; Maegio et al. 2000; Giidcl et afl 
12001; .Sanz-Forcada et al.. .2003; .Still fc Hussain ■20031) 
togeth er with our own. The values from IGiidel et alJ 
ll200l are from the RGS results. Note that the values 
for S and Si are few times smaller than those obtained 
with EPIC from the same observations. The S value 
obtained with EPIC is comparable with our values and 
previous works on AB Dor. The results obtained from 
the literature are all compared in Fig. ^ It is important 
to remember that different authors adopt different val- 
ues for "solar" photospheric abundances: for both the 
table and figure we have scaled the various abundance 
results to the photo spheric scale adopted here , based on 
the solar mixture of iGrevesse fc Sauvall ()1998() . which is 
similar to the stellar photospheric values derived for a 
few elements of the stars. Note that we have used the 
new solar photospheric oxygen a bundance, 0/H = 8.69, 
derived bv lAllende Prieto et all 1)20011) . 

While there is good agreement in the general trend of 
abundances as a function of FIP, there are clearly some 
significant differences between the different studies. The 
large spread in results from the different studies appear 
to be a result of modelling uncertainties and possibly to 
improvements of the atomic database. This conclusion 
is reinforced by the much more favourable comparison 
of the results of the recent Chandra and XMM-Newton 
studies of AB Do r — th is work compared with t hat of 
IGiidel et alJ (12001 and ISanz-Forcada et alJ l|200l. 

'Audard et al.H200 3all.'H uenemoerder et al.l(|2003|) and 
Sanz-Forc ada et all (2003) have recently reported, for 
V824 Ara, AR LAc and AB Dor respectively, higher 
values for abundances of Al, Ca, and Na than for ele- 
ments with slightly higher FIP suggesting a schism in 
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the strictly inverse FIP e ffect scfieme. In particular, 
iSanz-Forcada et aP l|200,'^ ) obtained different values of 
Ca abundance for AB Dor, including a very high value of 
[Ca/H]=0.62 derived from th e high resolution RGS spec- 
trum. Indeed, Figure 12 of 'Sanz- Forcada et alJ 1)200.'^ 
suggests that the inverse FIP effect in AB Dor under- 
goes a palpable reversal toward the lowest FIP elements. 
A sim ilar effect was suggested by iHuenemoerder et alJ 
lj2003j) based on high resolution Chandra spectra of 
AR Lac. That there might be a difference in abundance 
behaviour among the low FIP elements has already been 
suggested on the basis of solar measurements, in which 
the lowest FIP elements, such as Na and Al, appear 
to be enhanced relative to elements such as Si and Mg 
l)Feldm an'1992^. 

Although we have not been able to obtain definitive 
measurements of lines of Na, Al and Ca in our spectra 
of AB Dor (none were detected with any degree of sig- 
nificance), we have estimated approximate upper limits 
to the abundances of these elements through compari- 
son of observed and synthetic spectra. Fig 1101 illustrates 
the observed AB Dor MEG spectrum and two synthetic 
spectr a: one computed using the abundance derived for 
Ca byi anz-Forcada et al. (2003) and the other one as- 
suming Al, Ca and Na have coronal abundances similar 
to Mg. The former produce extremely strong Ca XIX, 
Al XII and Na XI lines which are clearly not seen in the 
observed spectrum. The synthetic spectrum computed 
assuming relative coronal abundances equal to Mg for 
Al, Ca and Na instead provide good approximations to 
the observed spectrum. We conclude that the coronal 
abundances of Al, Ca and Na in AB Dor are similar to, 
or lower than, those of the other low FIP elements Mg, 
Si and Fe. 

The coronal abundances we have derived for V471 Tau 
are a factor of two or less lower (factor thre e for Ne) than 
those reported bv fStill & HussainI l|2003D . This might 
be due to artifacts of the analysis of the low resolution 
ASCA pulse height spectra, though other explanations 
are possible. There were three modest flare events dur- 
ing the AS CA observ ation t hat w ere not filtered from the 
analysis of lStill fc HussainI l|2003j) . It appears likely that 
in some flares on active stars the abundances of low FIP 
elements are enhanced toward their pristine photospheric 
values, a result tentatively interpreted as chromosphcric 
evapor ation (e.g. reviews of Audard ct al. 2003b; Draki 
l2008aj) . Flares, then, might temporarily raise the abun- 
dances of the dominant plasma at a given time, though 
the ASCA flares were rather small for this to present an 
immediately convincing explanation. 

Alternatively, it is possible that the abundance pat- 
terns simply vary with time, according to the charac- 
teristics of dominant and relatively quiescent active re- 
gions, and at the time of the earlier ASCA observation 
the abundances were simply higher. Circumstantial ev- 
idence against such an explanation can be found in the 
abundance analysis of the active RS CVn-type binary 
HR 1099, for which independent analyses of Chandra and 
XMM-Newton observations obtained at different times 
yielded very similar values for qui escent corona l abun- 
dances (Brinkman et al. 2001; Drak e et al1l200lD . More- 
over, ^rgie (2003b) has found similar abundances in the 
coronae of HR 1099 and its "coronal twin" the active 
late-type component of the Algol binary. 



That our results for AB Dor and V471 Tau — coronal 
twins in many respects — are so similar reinforces the con- 
clusion that coronal abundances are dependent largely 
on stellar parameters and do not vary substantially dur- 
ing times of relative coronal quiescence. It also points 
out that tidal interactions, which may have an influence 
on differential rotation and possibly on dynamo action, 
are not important for determinnig coronal structure (the 
magnetic interactions between the white dwarf and the 
K star components of V471 Tau does not seem to play an 
important role at least in the X-ray domain). However, 
despite the very similar coronal temperature structure 
and abundance patterns, the Lx of V471 Tau is about 
twice higher than that for AB Dor (see Table This 
difference in brightness might be within the stochastic 
variation range for active stars, but might also be hint- 
ing that tidal interaction could be contributing to the 
stellar coronae somehow. 

6. CONCLUSIONS 

AB Dor and V471 Tau represent a rare pair stars that 
have the same rotational period {Porb ^ 0"?5) and spec- 
tral type but quite different evolutionary histories and 
ages(~70Myr vs '^GOOMyr). As such, a comparison 
between their coronal properties could provide an illu- 
minating glimpse of any fundamental underlying differ- 
ences in their magnetic dynamos and activity. Based on 
an analysis of high resolution Chandra X-ray spectra of 
these stars we draw the following conclusions. 

1. Both the coronal temperature structure and ele- 
ment abundances are remarkably similar. Element 
abundances are characterized by an inverse-FIP 
type effect in which the abundances of high FIP 
elements are enhanced relative to those of low FIP 
elements. Both coronae are depleted in low FIP 
elements relative to their photospheres by a factor 
of - 3. 

2. The observed abundance anomalies are less ex- 
treme than those typically seen in very active stars, 
despite AB Dor and V471 Tau lying at the bound- 
ary between coronal saturation and supersatur a- 
tion regimes l|Pizzolato et al.ll2003triardinell2004D . 

3. There is no evidence for any difference in the abun- 
dances of the very low FIP (— 5-6 eV) elements Na, 
Al and Ca as compared to the low FIP (~ 7.5-8 eV) 
elements Si, Mg and Fe, as some earlier work has 
suggested. 

4. Abundances derived using our method in which the 
DEM is obtained from a combination of the ratios 
of a few well-selected prominent spectral lines and 
continuum measurements are robust and indepen- 
dent of both the exact set of lines used to define the 
DEM and of the details of its temperature struc- 
ture. 

5. Small-scale structure is found in the DEMs ob- 
tained based on different sets of spectral lines. Such 
structural details are most likely spurious artifacts 
arising from either errors in assumed ion popu- 
lations or spectral line excitation rates, or non- 
uniform abundances in different coronal structures 
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that contribute significantly to the observed X-ray 
emission. Coronal DEMs are strongly dependent on 
the input data used to derive th em, as was empha- 
sised e arlier in the solar study of lKashvap fc Drak3 

6. Based on the very similar coronal characteristics of 
AB Dor and V471 Tau we conclude that there is no 
observable difference between the coronal activity 
of youth and high primordial angular momentum 
and that of tidal spin-up of an older star. In other 
words, the exact evolutionary state has little effect 
on the corona, and the parameters that dominate 
coronal structure and composition are rotation rate 
and spectral type. The latter determines of course 
the parameters of the convection zone, where the 
magnetic fields sustaining coronal activity are be- 
lieved to originate, while rotation rate appears to 
be the dominant parameter in the amplification of 
these magnetic fields. 

7. We also note that despite the binarity of V47I 
Tau it does not show any significant difference in 
thermal structure and composition with respect to 
AB Dor. This would mean that tidal interactions, 
which may have an influence on differential rotation 
and possibly on dynamo action, are not important 
for determining coronal structure. The magnetic 
interactions between the white dwarf and the K 
star components of V471 Tau does not seem to 
play an important role at least in the X-ray do- 
main. The latter conclusion was previously sug- 
gested in an ana lysis based on radio observations 
l|Lim et aLlll99(iD . Nevertheless, the observed dif- 
ference in brightness between AB Dor and V471 
Tau (twice higher) might be hinting that tidal in- 
teraction could be contributing to the stellar coro- 
nae somehow. 
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TABLE 1 

Summary of Stellar Parameters. 





Spectral 


Period 


Distance 


Active 


Component 




V sini 




Star 


Type 


[d] 


[pc] 


Te//[K] 


Re 


[deg] 


[kms-i] 


[1030 orgs-1] 


AB Dor 


KO-2 V 


0.51*^ 


14.9'^ 


5250'* 


0.86° 0.76^= 


eo' 


93'' 


0.70' 


V471 Tau 


K2V + WDS 


0.52'^ 


48' 


5040j 


0.96j 0.93j 


77j 


91"^ 


1.58' 



^Parameters of the component that d ominates X-ray activity: AB Dor and the K2 V component of V471 Tau. 
'^Collier Ca meron &: Robinsonl 1198911 
tfrakull I 198^) 
"[Guirado et^Q 119971) 

"Johnson 1 196^) 

Masigio ct al. 



^Cucrster ct al. 
{Rucinski 1 1981 


1 I99I1 

) 


'^uinan & Ribas 




|de Bruiine ct al. 


(2001) 



adopted; in particular the mass and radius arc larger for smaller inclinations. 
" iRamsever et al.l fl995) 
'Derived from our Chandra observations. 



TABLE 2 
Summary Chandra Observations. 





ObsID 


Instrument 


Start 


Exp 


Star 






[UT] 


[ks] 


AB Dor 


0016 


HETC-I-ACIS-S 


1999-10-09T11:21:09 


52.3 


V471 Tau 


2523 


LETC'^-I-HRC-S 


2002-01-24T22:30:03 


87.5 



IWeisskopf et al.l i20o3) 



' iPredehl et al.l 1199711 
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TABLE 3 

Identification and fluxes for spectral lines, observed on AB Dor and V471 Tau, used in this analysis. 



'^obs '^pred 

(A) (A) 



Ion 



(K) 



AB Dor 
(10-5 ph/cm2/s) 



V471 Tau 
(10-5 ph/cm2/s) 



Transition 
upper — » lower 



4.722 


4.727 


S XVI 


7.40 


1.56 ±0.59 


0.33 ± 0.72 




4.733 


S XVI 


7.40 






5.032 


5.039 


S XV 


7.20 


2.26 ± 0.75 


0.62 ± 0.70 


6.177 


6.180 


Si XIV 


7.20 


4.80 ± 0.38 


0.59 ±0.28 




6.186 


Si XIV 


7.20 






6.647 


6.648 


Si XIII 


7.00 


4.79 ± 0.40 


0.85 ± 0.34 


8.422 


8.425 


Mg XII 


7.00 


5.95 ± 0.39 


0.88 ± 0.35 




8.419 


Mg XII 


7.00 






9.173 


9.169 


Mg XI 


6.80 


3.71 ±0.39 


0.61 ±0.42 


10.028 


10.029 


Na XI 


6.90 


2.89 ± 0.65 


0.73 ± 0.30 




10.023 


Na XI 


6.90 






10.998 


11.003 


Na X 


6.70 


8.56 ± 0.85 


1.23 ±0.36 


12.133 


12.132 


Ne X 


6.80 


68.06 ± 1.80 


8.00 ± 0.77 




12.137 


Ne X 


6.80 






12.283 


12.285 


Fe XXI 


7.00 


14.33 ± 1.33 


1.74 ±0.52 


13.433 


13.447 


Ne IX 


6.60 


32.96 ± 1.79 


3.02 ± 0.53 


13.508 


13.504 


Fe XIX 


6.90 


13.87 ± 1.62 


3.58 ±0.57 


14.208 


14.208 


Fe XVIII 


6.90 


21.24 ± 1.62 


2.29 ±0.54 




14.203 


Fe XVIII 


6.90 






14.263 


14.267 


Fe XVIII 


6.90 


5.69 ± 1.09 


1.27 ±0.47 


15.013 


15.015 


Fe XVII 


6.75 


36.18 ±2.01 


5.23 ± 0.62 


16.008 


16.007 


O VIII 


6.50 


33.87 ±2.13 


3.56 ± 0.60 




16.006 


O VIII 


6.50 






18.973 


18.973 


O VIII 


6.50 


154.21 ±6.18 


14.32 ± 1.05 




18.967 


O VIII 


6.50 






21.607 


21.602 


O VII 


6.30 


17.49 ±4.16 


2.71 ± 0.84 


24.777 


24.779 


N VII 


6.30 


23.01 ±4.70 


4.84 ± 1.07 




24.785 


N VII 


6.30 







(2p) 2P3/2 ^ 
(2p) 2P1/2 ^ 
(ls.2p) IPi ^ 
(2p) 2P3/2 -» 
(2p) 2P1/2 ^ 
(ls.2p) IPi ^ 
(2p) 2P1/2 ^ 
(2p) 2P3/2 ^ 
(ls.2p) IPi 
(2p) 2P1/2 ^ 
(2p) 2P3/2 ^ 
(ls.2p) ll'i ^ 
(2p) 2P3/2 ^ 
(2p) 2P1/2 -> 
(2s2.2p3d) 3Di 



(Is) 2S1/2 
(Is) 2S1/2 
(ls2) ISo 
(Is) 2S1/2 
(Is) 2S1/2 
(ls2) ISo 
(Is) 2S1/2 
(Is) 2S1/2 
(ls2) iSo 
(Is) 2S1/2 
(Is) 2S1/2 
(ls2) ISo 
(Is) 2S1/2 
(Is) 2S1/2 

(2s2.2p2) 3Po 



(ls.2p) IPi (ls2) ISo 
(2p3(2P).3d) ID2 (2s2.2p4) 3P2 
(2p4(lD).3d) 2P3/2 ^ (2s2.2p5) 2P3/2 
(2p4(lD).3d) 2D5/2 ^ (2s2.2p5) 2P3/2 
(2p4(lD).3d) 2Fs/2 ^ (2s2.2p5) 2P3/2 
(2p5.3d) IPi (2p6) ISo 



(3p) 2P1/2 - 
(3p) 2P3/2 - 
(2p) 2P1/2 - 
(2p) 2P3/2 - 
(ls.2p) IPi 
(2p) 2P3/2 - 
(2p) 2P1/2 - 



(Is) 2S1/2 

(Is) 2S1/2 
(Is) 2S1/2 
(Is) 2S1/2 
* (ls2) ISo 
(Is) 2S1/2 
(Is) 2S1/2 



^The temperature at which the function Gij(T) peaks. 
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TABLE 4 

Coronal abundances obtained from an abundance-independent DEM-reconktruction for AB Dor and V471 Tau, relative to stellar 

PHOTOSPHERIC VALUES. NOTE: THE ABUNDANCES HAVE BEEN CALCULATED FROM FOUR DIFFERENT SET OF IONS FOR EACH OBJECT. 



Element" FIP'' X/H" AB Dor'' AB Dor" AB Dor' AB Dor" 



[Na/H] 


5.14 


6.33 


[Al/H] 


5.98 


6.47 


[Ca/H] 


6.11 


6.36 


[Mg/H] 


7.65 


7.58 


[Fe/H] 


7.87 


7.50 


[Si/H] 


8.15 


7.55 


[S/H] 


10.36 


7.33 


[O/H] 


13.62 


8.69 


[Ne/H] 


21.56 


8.08 



<0.54'' <0.54'' 
<0.54'' <0.54'' 
<0.54'' <0.54'' 
-0.56± 0.03 -0.56± 0.03 
-0.57± 0.04 -0.56± 0.04 
-0.50± 0.03 -0.49± 0.04 
-0.51± 0.11 -0.49± 0.11 
-0.22± 0.04 -0.21± 0.04 
-0,10 ' 0.03 -0,09-^ 0.03 



<0.54'' <0.5l'' 

<0.54'' <0.5l'' 

<0.54'' <0.5l'' 

-0.61± 0.03 -0.60± 0.03 

-0.47± 0.03 -0.54± 0.03 

-0.53± 0.03 -0.50± 0.03 

-0.54± 0.11 -0.47± 0.11 

+0.07± 0.04 -0.14± 0.04 

-0.03 f 0.03 -0.06+ 0.03 



Element'^ FIP'" X/H° V471 Tau'^ V471 Tau° V471 Tau'' V471 Tau^ 



[Na/H] 


5.14 


6.33 


[Al/H] 


5.98 


6.47 


[Ca/H] 


6.11 


6.36 


[Mg/H] 


7.65 


7.58 


[Fe/H] 


7.87 


7.50 


[Si/H] 


8.15 


7.55 


[S/H] 


10.36 


7.33 


[O/H] 


13.62 


8.69 


[Ne/H] 


21.56 


8.08 



-0.58± 0.17 -0.56± 0.17 
-0.70± 0.12 -0.66± 0.11 
-0.53± 0.15 -0.48± 0.15 
-0.37± 0.46 -O.Slib 0.45 
-0.32± 0.10 -0.32± 0.10 
-0.24± 0.09 -0.24± 0.09 



-0.57± 0.16 -0.58± 0.16 
-0.64± 0.09 -0.66± 0.09 
-0.48± 0.15 -0.50± 0.15 
-0.30± 0.45 -0.33± 0.45 
-0.25± 0.10 -0.29± 0.09 
-0.21± 0.09 -0.23± 0.08 



"Abundances relative to stellar photosplieric values. 
'First Ionization Potential in eV. 

'^Assumed pliotospheric abundances expressed relative to H by number on a logarithm scale where H/H=12. 

''Abundances based on DElVIs from O, Ne, Mg and Si. 
'^Ditto from O, Ne, Mg, Si and S. 
^Ditto from O, Si, S, FeXVII, FeXVIII and FeXXI. 
"Ditto from O, Ne, Mg, Si, FeXVII, FeXVIII and FeXXI. 

''Upper limits based on synthetic spectra estimation. 



TABLE 5 

AB Dor and V471 Tau Abundance Ratios using Temperature-Insensitive Diagnostics. 



Abundance Ratio AB Dor V471 Tau 



[N/O] 0.25± 0.09 0.53± 0.11 

[O/Ne] -0.23± 0.11 -0.19± 0.15 

[Ne/Mg] 0.47± 0.04 0.32± 0.24 

[Nc/Fe] 0.71± 0.03 0.57± 0.07 

[Mg/Si] -0.09± 0.04 -0.09± 0.22 

[Mg/Fe] -0.18± 0.05 0.03± 0.31 

[Si/S] O.Olifc 0.12 -O.lSifc 0.47 



The Coronae of AB Dor and V471 Tau 



15 



TABLE 6 

Comparison of coronal abundances from AB Dor and V471 Tau obtained from this work and the literature. Note that all the coronal 

ABUNDANCES ARE RELATIVE TO STELLAR PHOTOSPHERIC VALUES. 



Element^ 


Fipb X/H<: AB Dor'* 


AB Dor<= 


AB Dor' 


AB DorS 


AB Dorh 


AB Dor* 


V471 Tauj 


V471 Tau* 


[Na/H] 


5.14 


6.33 












-0.58"^ 






[Al/H] 


5.98 


6.47 








-0.62 


-0.35 


-0.58'' 






[Ca/H] 


6.11 


6.36 


-0.31 




-0.74 


-0.99 


0.41 


-0.58"^ 






[Ni/H] 


7.63 


6.25 


-0.35 


-0.36 


-0.33 


-0.64 


-0.18 








[Mg'/H] 


7.64 


7.58 


-0.24 


-0.72 


-0.57 


-0.58 


-0.39 


-0.58 


-0.06 


-0.57 


[Fc/H] 


7.87 


7.50 


-0.47 


-0.73 


-0.48 


-0.59 


-0.40 


-0.54 


-0.02 


-0.67 


[Si/H] 


8.15 


7.55 


-0.46 


-0.80 


-0.85 


-0.50 


-0.41 


-0.51 


-0.12 


-0.50 


[S/H] 


10.36 


7.33 


-0.55 


-0.38 


-1.52 


-0.66 


-0.46 


-0.50 




-0.33 


[C/H] 


11.26 


8.52 






-0.30 


-0.32 


-1-0.00 








[O/H] 


13.62 


8.69 


-0.14 


-0.26 


-0.16 


-0.31 


-1-0.16 


-0.13 


-1-0.19 


-0.30 


[N/H] 


14.53 


7.92 


-0.32 




-0.15 


-0.21 


-1-0.22 








[Ar/H] 


15.76 


6.40 


-1-0.16 




-f0.09 


-0.30 


-f0.19 








[Ne/H] 


21.56 


8.08 


-1-0.09 


-0.20 


+0.01 


-0.11 


-t-0.23 


-0.07 


0.39 


-0.23 



"Abundances relative to stellar photospheric values. 
*First Ionization Potential in eV. 

'^Assumed photospheric abundances expressed relative to H by number on a logarithm scale where H/H=12. 
Mewe et ail iW9^ 
Ortolani et al.l il9M) 

Eiidel et al.' ("2OOI') based on RGS resuhs 
mz-Forcada et al. (2003,) based on 3-T analysis. 
1Sanz-Forcada et alJ 120031) based on DEM analysis. 

'itiS-Wflrk , , . 

^StiU & H ussaiiJ IpOOS) 

''Upper limits based on synthetic spectra estimation. 
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-0 VI 
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Olo o> 00 o to Olo 01 00 o to 
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Fig. 1. — Chandra X-ray spectra of AB Dor (black shade lower panel; HETG-I-ACIS-S) and V471 Tau (grey shade upper panel; 
LETG-I-HRC-S) respectively. The strongest lines over the observed wavelength range are identified. The thick black line in the upper 
panels shows AB Dor spectra after smoothing with a beta-profile in order to match with the resolution of the V471 Tau spectrum. Note 
the strong correspondence between the two spectra. 
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Fig. 2.— Chandra X-ray light curves of AB Dor (top; HETG+ACIS-S) and V471 Tau (bottom; LETG-I-HRC-S) respectively. The AB 
Dor and V471 Tau light curves are bin ned at 100s and 200s int ervals respectively. Both objects were relatively quiescent, showing no large 
flare events. Based on the ephemeris of lGuinan Hz Ribad |20^, the V471 Tau light curve spans phases (p = 0.05 to 1.97. 
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Fig. 3. — Comparison of observed and blend corrected fluxes of AB Dor and V471 Tau Chandra X-ray spectra. Note that the fluxes are 
corrected for distance. 
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Fig. 6. — Comparison between LETG+HRC-S spectrum of V471 Tau (top panel) and the predicted spectrum derived from the four 
reconstructed DEMs. 
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Fig. 7. — Top panels: AB Dor DEM obtained by running a MCMC[M] reconstruction code on a set of lines of H-likc, He-like and highly 
ionized Fe line fluxes. The panels from left to right correspond to line set 1 (O, Ne, Mg, Si), set 2 (O, Nc, Mg, Si and S), set 3 (O, Si, S, 
Fe XVII, Fc XVIII and Fe XXI) and set 4 (O, Ne, Mg, Si, Fc XVII, Fc XVIII and Fc XXI). The thick solid line represents the best-fit 
DEM, while the shaded regions correspond to the l-cr deviations present in each temperature bin. Bottom panels: Same as top panels but 
for V471 Tau. 
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Fig. 9. — Comparison of coronal abundances against FIP from AB Dor and V471 Tau obtained from this work (black symbols) 
and the lite rature (grey s ymbo ls). The uncertainties in the abundances derived in this work are of order 0.1 dex. Note that the 
values from lGiidel et ahl (|2001f) are from the RGS results. 
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